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Abstract 

'-^ ■ Vertically coupled double quantum rings submitted to a perpendicular magnetic field B are 

addressed within the local spin-density functional theory. We describe the structure of quantum 



ring molecules containing up to 40 electrons considering different inter-ring distances and intensities 
of the applied magnetic field. When the rings are quantum mechanically strongly coupled, only 



bonding states are occupied and the addition spectrum of the artificial molecules resembles that 

a, 

of a single quantum ring, with some small differences appearing as an effect of the magnetic 



field. Despite the latter has the tendency to flatten the spectra, in the strong coupling limit some 



■ clear peaks are still found even when B ^ that can be interpretated from the single-particle 

m ; 

energy levels analogously as at zero applied field, namely in terms of closed-shell and Hund's-rule 
configurations. Increasing the inter-ring distance, the occupation of the first antibonding orbitals 
washes out such structures and the addition spectra become flatter and irregular. In the weak 

00 ' 

O ■ coupling regime, numerous isospin oscillations are found as a function of B. 

> '■ 

^ . PACS numbers: 85.35.Be, 73.21.-b, 71.15.Mb, 75.75.+a 
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I. INTRODUCTION 



Systems made of correlated electrons confined in semiconductor nanoscopic dot and ring 
structures, so-called quantum dots (QDs) and rings (QRs) respectively, have been the subject 
of intense theoretical and experimental research, see e.g. Refs. |l|j2| and references therein. 
From the latter point of view, for quantum dots it has been proved^ the possibility to tune 
over a wide range the number of electrons contained in the system, as well as to control 
both the size and the shape of the dots by means of external gate voltages, goal that has 
not been achieved yet for ring geometries due to the higher complexity of their fabrication 
process,-^ which involves several experimental techniques such as atomic force microscopy- 
strain-induced self-organization^ or droplet molecular beam epitaxy. 8 

The interest of QRs arises from their peculiar behavior in the presence of a perpendicularly 
applied magnetic field (B), which is very distinct from that observed in QDs and shows 
up as an oscillatory behavior of their energy levels as a function of B. This property, 
together with the fact that in narrow enough QRs the electrons experiment a nearly one- 
dimensional Coulomb repulsion, leads to the integer and fractional Aharonov-Bohm effects, 
usually associated with the appearance of the so-called persistent currents in the ring.- These 
quantum-interference phenomena have been experimentally reported 10 and have motivated 



a series of theoretical works whose number is steadily increasing, see e.g., Refs. 
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16l and references therein. 
One of the most appealing possibilities offered by electron systems confined in semi- 
conductor heterostructures is their ability to form coupled entities, usually referred to as 
"artificial molecules" , in which the role of the constituent "atoms" is played by single quan- 
tum dots or rings and that have analogies with natural molecules such as the hybridization 
of the electronic states forming molecular-like orbitals. In addition, these artificially cou- 
pled systems present important advantages such as a tunable "interatomic" coupling by 
means, e.g., of the modification of the relative position/size of the constituents. This fact 
has, besides its intrinsic interest, potential relevance to quantum information processing 
schemes since basic quantum gate operations require controllable coupling between qubits. 
In this sense, artificial molecules based on two coupled QDs, called quantum dot molecules 
(QDMs) have been proposed as scalable implementations for quantum computation pur- 
poses and have received great attention from the scientific community in the last years -see 
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e.g. Refs. 0QQ|2O|j2l|j2^j23|j24j25 



|26| and references therein. 
Also, molecular-beam epitaxy techniques have recently allowed the synthesis of quantum 
ring molecules (QRMs) in the form of concentric double QRs^>2£ and vertically stacked 
layers of self-assembled QRs,— ^ the optical and structural properties of the latters having 
also been characterized by photoluminescence spectroscopy and by atomic force microscopy, 
respectively. This has sparked theoretical studies on the structure and optical response of 
both vertically and concentrically coupled QRs of different complexity and scope, revealing 
properties different from those of their dot counterparts due to the non-simply connected 
ring topology. For instance, studies on the single-electron spectrum of vertical QRMs^ 1 ^ 
have shown that the electronic structure of these systems is more sensitive to the inter-ring 
distance than that of coupled QDs. As a consequence, in ring molecules quantum tunneling 
effects are enhanced since less tunneling energy is required to enter the molecular-type 
phase. Also, the consideration of "heteronuclear" artificial molecules constituted by slightly 
different QRs offers the interesting possibility to control the effective coupling of direct- 
indirect excitons^ by means of the application of a magnetic field and taking advantage of 
the fact that charge tunneling between states with distinct angular momentum is strongly 
suppressed by orbital selection rules. To this end, some authors have considered the case of 
QRMs made of strictly one- dimensional, zero-thickness QRs and have used diagonalization 
techniques to address the few-electron problem.— £iL2Mk The simultaneous effect of both 
electric and magnetic fields applied to a single-electron QRM has also been studied^ -see 



also Ref. 
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and the optical response of QRMs where the thickness of the constituent 
QRs is taken into account has been obtained.— In addition, the spatial correlation between 
electron pairs in vertically stacked QRs only electrostatically coupled has been shown to 
undergo oscillations as a function of the magnetic flux, with strongly correlated situations 
between ground states with odd angular momentum turning out to occur even at large 
inter-ring distances.— More recently, the structure of a QRM made of two vertically stacked 
quantum rings has been addressed at zero magnetic field for a few tens of electrons within the 
local spin-density functional theory (LSDFT) neglecting^ and incorporating^ the vertical 
thickness of the constituent QRs. 

In this work we address the ground state (gs) of two thick, vertically coupled identical 
quantum rings forming "homonuclear" QRMs populated with up to 40 electrons and pierced 
by a perpendicularly applied magnetic field. We extend in this way our previous study— 
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addressing the appearance and physical interplay between the spin and isospin^ 3 . degrees 
of freedom as a function of the variation of both the intensity of the magnetic field and 
the inter-ring separation. Modelling systems charged with such large number of electrons 
requires the employment of methodologies that minimize the computational cost. Here we 
have made use of the LSDFT,— ^ whose accuracy for the considered values of the magnetic 
field has been assessed 24 by comparing the obtained results for a single QD with those given 
by the current-spin-density functional theory (CSDFT),— which in principle is better-suited 
for high magnetic fields, and also with exact results for artificial molecules.— 

This paper is organized as follows. In Sec. II we briefly introduce the LSDFT and the 
model used to represent the vertical QRMs. In Sec. Ill we discuss the obtained results for 
some selected configurations, and a summary is given in Sec. IV. 



II. DENSITY FUNCTIONAL CALCULATION FOR MANY-ELECTRON VERTI- 
CAL QUANTUM RING HOMONUCLEAR MOLECULES 

The axial symmetry of the system allows one to work in cylindrical coordinates. The 
confining potential V c f(r,z) has been taken parabolic in the xy-pl&ne with a repulsive core 
around the origin, plus a symmetric double quantum well in the ^-direction, each one with 
width w, depth Vo, and separated by a distance d. To improve on the convergence of the 
calculations, the double-well profile has been slightly rounded off, as illustrated in Fig. 2 of 
Ref. I24J. The potential thus reads V c f(r,z) = V r (r) + V z (z), where 

V r (r) = V Q(R -r)+ l -mul{r- R ) 2 0(r - R ) 

VU) = V < ±+e {z+ * /2+w)/ ' T ~ l+ e (2+<V2)/<7 if ^ < 

[ l +e (z-d/2)/ CT ~ 1+e (z-d/2-w)/o- ^ Z > , 

with cr = 2 x 10~ 3 nm, and Q(x) = 1 if x > and zero otherwise. The convenience of using 
a hard-wall confining potential to describe the effect of the inner core in QRs is endorsed 
by several works in the literature.— We have taken Ro = 10 nm, Vo = 350 meV, Tiujo = 6 
meV and w = 5 nm. These parameters determine the confinement for the electrons together 
with the distance between the constituent quantum wells that is varied to study QRMs in 
different inter-ring coupling regimes. 

For small electron numbers (iV), it is justified to take u)q to be iV-independent. However, 
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in a more realistic scheme its value should be tuned according to the number of electrons 
contained in the system, relaxing the confinement as the latter is increased. In the case of 
quantum dots it has oftenly been used a iV _1//4 -dependence that arises from the r-expansion 
near the origin of the Coulomb potential created by a two-dimensional uniform positive 
charge distribution -jellium model- and that it is generalized to the case of quantum dot 
molecules as u>o = ftiV B 1//4 , Nb being the number of electrons filling bonding orbitals -see 
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It is clear that the mentioned 



below. The rationale for this generalization is given in Ref. 
iV-dependence would be harder to justify for QRs, and in fact no alternative law is known 
for a single QR that could be generalized to the case of QRMs. For this reason, in this work 
we have taken u to be iV-independent, which is to some extend less realistic for the largest 
values of N we have considered. 

Considering the iV-electron system placed in a magnetic field parallel to the z-axis, 
within the LSDFT in the effective mass, dielectric constant approximation, the Kohn-Sham 
equations^ 4 - 1 ^ in cylindrical coordinates read 

I ( 8 2 Id l 2 d 2 \ U Cl 1 2 2 Tr . . 



+ Vu+ r„ + ( W xc + -g*fJL B B) >ia 



(2) 

U n i a {r,z) = £ n laU n l a {r,z) , 



where the single-particle (sp) wave functions have been taken to be of the form 
(finiair, z, 8, a) = u n i a (r, z)e~ tW Xa with n = 0, 1, 2, . . ., I = 0, ±1, ±2, . . ., -I being the pro- 
jection of the single-particle orbital angular momentum on the symmetry axis, and er=T(!) 
representing spin-up(down) states. The vector potential has been chosen in the symmetric 
gauge, namely A = B{— y, x, 0)/2; \ib = he/(2m e c) and u c = eB/c are, respectively, the 
Bohr magneton and the cyclotron frequency, and i] a =+l(— 1) for cr=|(|); Vn(r,z) is the 
direct Coulomb potential, and V xc = d£ xc (n,m)/dn\ gs and W xc = dS xc (n,m)/dm\ gs are 
the variations of the exchange-correlation energy density £ xc (n,m) in terms of the electron 
density n(r,z) and of the local spin-magnetization m(r,z) = n^(r,z) — n^(r,z) taken at 
the gs. £ xc (n,m) = £ x (n,m) + £ c (n,m) has been built from three-dimensional homoge- 
neous electron gas calculations; this yields a well-known,— simple analytical expression for 
the exchange contribution £ x (n,m). For the correlation term £ c (n,m) we have used the 
parametrization proposed by Perdew and Zunger™< Details about how the Kohn-Sham and 



the Poisson equations have been solved can be found in Ref. 
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Notice the use in Eq. (j2J) 



of effective atomic units H = e 2 /e = m =1, where e is the dielectric constant and m the 
electron effective mass. In units of the bare electron mass m e one has m = m*m e , the length 
unit being the effective Bohr radius = aoe/m* and the energy unit the effective Hartree 
H* = Hm*/e 2 . In the numerical applications we have considered GaAs quantum rings, for 
which we have taken e = 12.4, and m* = 0.067; this yields a* Q ~ 97.9 A and H* ~ 11.9 meV, 
the effective gyromagnetic constant being g* = —0.44. 

To label the gs configurations ("phases") we use an adapted version of the ordinary 
spectroscopy notation,— namely 2S+1 L^ U , where S and L are the total and \L Z \, re- 
spectively. The superscript +(— ) corresponds to symmetric (antisymmetric) states under 
reflection with respect to the z = plane bisecting the QRMs, and the subscript g(u) refers 
to positive (negative) parity states. All these are good quantum numbers even in the presence 
of an axial magnetic field. By analogy with natural molecules, symmetric and antisymmet- 
ric states are referred to as bonding (B) and antibonding (AB) orbitals, respectively. We 
have defined the "isospin" quantum number I z -bond order in Molecular Physics- as^^ 1 ^ 
I z = (Nb — Nab)/2, Nb(ab) being the number of occupied bonding(antibonding) sp states. 

III. RESULTS 

Due to the large number of variables needed to characterize a given QRM configuration 
(electron number, magnetic field and inter- ring distance), we limit ourselves to present results 
in a limited range of values for such variables, aiming at discussing calculations that might 
illustrate the appearance of some properties of the systems under study. For the sake of 
comparison, we have also addressed one single QR symmetrically located with respect to 
the z = plane with the same thickness (5 nm) and radial confinement as the coupled rings. 

Fig. [TJ shows the Kohn-Sham sp levels for one single ring hosting A=40 electrons as a 
function of I for different values of the applied magnetic field. As it is well known, these 
levels are ±/-degenerate at B — 0. In this particular case, the gs has S z = 1, and it is made 
up of symmetric (with respect to z = 0) sp states with up to n = 3. In the non-interacting 
single-electron model, in which the Coulomb energy is not considered and consequently 
the sp wave functions factorize into a r-dependent and a ^-dependent part with associated 
quantum numbers n r and n z , i.e., u n i(r, z) — > U nr (r)Z nz (z), one would say that the gs is 
made up of sp states with n z = and radial quantum numbers up to n r = 3. 
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When B ^ 0, the ±/- degeneracy is lifted and, on the other hand, the I < sp levels 
become progressively depopulated in favor of those with I > as the magnetic field increases 
until eventually -at about ~ 4 T- only I > orbitals are filled. At this point, only a few 
states with n = 2 are occupied, and the ring has S z = 0. From this value of B on, the 
simultaneous filling of increasingly higher-/ states and those close to / = gives rise to 
configurations containing only states with n = 1 and with large values of the total spin 
(e.g., S z = 9 for B = 8 T). Eventually, the system becomes fully spin-polarized at B ~ 13.5 
T. It is worth noticing the conspicuous bending of the "Landau bands" (sets of bonding or 
antibonding states characterized by the same n and spin, and different value of /), instead of 
displaying a fairly flat region, as it happens when the in-plane confinement is produced by 
a jellium-like potential,— but not with our present choice of a A-independent parabola. It 
is also worth to stress that, due to the much stronger confinement in the vertical direction 
as compared to that in the radial one, only symmetric states are occupied. 

Analogously, the energy levels corresponding to QRMs with N = 40 and inter-ring dis- 
tances d — 2, 4 and 6 nm are shown in Figs. EH One can see the gradual evolution of the 
system as d increases; indeed, at d = 2 nm the spectrum is very similar to that of the single 
ring, with only bonding sp states being occupied. As d increases, a few antibonding orbitals 
become populated at small 5's, as one can see from the top panels of Fig. [3j corresponding 
to d — 4 nm, but eventually for increasing values of B the QRMs have again ground states 
where only bonding states are populated, as can be seen from the bottom panels of the same 
figure. For this inter- ring distance, the fully spin-polarized state is reached at B ~ 13.75 T. 
Finally, for the largest ring separation considered, namely d = 6 nm, a large amount of anti- 
bonding orbitals become occupied giving rise to small J z 's instead of the fairly large isospin 
values found for similar configurations at smaller distances (compare the bottom panels in 
Fig. |4] with those in Figs. [2] and [3]). In particular, the fully spin-polarized gs is found at 
about B ~ 7 T with I z = 2, whereas for d = 2 and 4 nm it appears near B = 14 T and has 
the maximum possible isospin value, namely I z = 20. At d — 6 nm, the maximum-spin state 
naturally consists of two distinct bands, one made up of bonding and another of antibond- 
ing states. These configurations are the QRM-analogues of the maximum density droplet 
(MDD) configurations found for QDMs at similar inter-dot distances, called respectively 



MDD B and MDD AB in Ref. 
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Increasing further the magnetic field causes the progressive 



occupation of higher-/ orbitals, which provokes the depopulation of the antibonding band 
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and the consequent increase of I z . For the highest considered magnetic field (B ~ 14 T), 
some antibonding orbitals are still occupied, yielding I z — 17. 

These results are a consequence of the evolution with d of the energy difference between 
bonding and antibonding states, AgASi which accurately varies as a function of the inter- 
ring distance according to the law As as = A e~ d ^ d °, already found for QDMs. 22 In our case, 
from the difference in energy of single-electron (bonding and antibonding) QRMs we have 
obtained Ao = 82 meV and do = 1.68 nm,— values which have turned out to be unaffected by 
the applied magnetic field. Clearly, the value of hu as compared to A S as, which allows to 
discern between the strong (hu < As as) and the weak (hu ^> A^s) quantum mechanical 
coupling regimes, has a crucial influence on the actual filling of bonding and antibonding sp 
states at a given inter-ring distance. Indeed, increasing ou while keeping constant the double 
well structure may favor the population of antibonding orbitals for large enough values of 
N .— This can be understood from the non-interacting electron model, in which the single- 
electron energies are the sum of two independent terms, one arising from the ^-localization 
and characterized by the quantum number n z , and another, which increases as ujq does, 
arising from the r-localization and depending on / and the radial quantum number n r . UN 
is large enough, the QRMs can minimize its energy by populating antibonding states with 
low values of n r and I instead of going on populating bonding states with higher quantum 
numbers. This explains why some antibonding states were filled even for d = 2 nm at B = 



and N = 40 for the QRMs of Ref. |39|, where ujq was taken to be 15 meV, value almost three 
times larger than the one considered in the present work. 24 

This particular structure of the bonding and antibonding bands at high magnetic fields 
may have some observable effects on the far-infrared response of QRMs. Indeed, since the 
dipole operator cannot connect bonding with antibonding sp states, for QRMs in the weak 
coupling limit one would expect the dipole spectrum to display additional fragmentation 
in the characteristic edge modes of the ring geometry^ due to the contribution of the 
antibonding electron- hole pairs (see e.g. the bottom panels of Fig. Hj). 

Figure [5] shows the evolution with d of the gs energy and the molecular phase of QRMs 
made up of N = 8 electrons and submitted to magnetic fields of different intensities. Notice 
that even moderate values of B give rise to ground states with large total angular momentum, 
which increases as the magnetic field does. For this reason, we have denoted it by its actual 
value instead of employing the usual notation with upper Greek letters, except for the cases 
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with L z = 0. Similar conclusions can be drawn for all the values of the magnetic field we 
have considered; on the one hand, for the studied inter-ring distances, the energy of the 
molecular phases increases with d due to the enhancement of the energy of the bonding 
states,— which dominates over the decrease of the Coulomb energy -for larger distances the 
constituent QRs are so apart that eventually this decrease dominates and the tendency is 
reversed. On the other hand, one can see that the first phase transitions are always found at 
the largest inter-ring distances since, as happens for QDMs in the few-electron limit, they 
are due to the replacement of an occupied bonding sp state by an empty antibonding one. 
This also explains why in most of the cases, and especially for the highest magnetic fields, 
the total angular momentum of the QRMs in the weakest coupling regime is reduced: the 
filled antibonding orbitals have lower Z's than the replaced bonding states. 

We have determined the magnetic field that gives rise to ring molecules with fully spin- 
polarized gs, and show it in Fig. [6] as a function of TV for different inter-ring distances going 
from the strong to the weak quantum mechanical coupling regimes. The isospin value of 
each configuration is also indicated. The number of electrons, N = 8 x M with M — 1 
to 5, was chosen with the aim of producing closed-shell structures at B = in the weak 
coupling limit. One can see that the results for d = 2 and 4 nm are very close, with only 
noticeable differences for N = 32. This can be understood from the bottom panels of Figs. 
|2] and [31 which show that for rather large magnetic fields only bonding orbitals are occupied 
for both ring separations. Contrarily, from Fig. H] one can see that in weaker coupling 
regimes the filling of antibonding states favors the fully spin-polarization of the QRMs at 
low B intensities as compared to those needed when the rings are closer to each other, which 
explains the differentiated results corresponding to d = 5 and 6 nm in Fig. [61 

When antibonding orbitas are populated, the variation of the magnetic field yields nu- 
merous transitions between different molecular phases with different isospin that are more 
complex than these observed in vertically coupled QDs. This particular behavior is mainly 
due to the periodic destabilization suffered by the lowest-Z occupied orbitals induced by the 
magnetic field, which is a direct consequence of the Aharonov-Bohm effect and makes it 
rather difficult to find a pattern among the observed evolutions for the different electronic 
populations. The spin and isospin phases as a function of the magnetic field are shown in 
Fig. [7] for d = 6 nm corresponding to N = 8, 16, and 24. It can be seen that in all cases 
at B = the QRMs have I z — 2 and S z — 1; when B is increased, non-monotonic spin 
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and isospin oscillations with AI Z = ±1 and AS Z = ±1 and 2 appear, respectively. Two 
facts, also present in QDMs,— ^ are worth to be stressed: on the one hand, molecular phase 
changes from — (+) to +(— ) ground states -recall that, as explained in Sec. II, this sign 
is related to the symmetry of the molecular configuration- involve AI Z = +1(-1) flips; on 
the other hand, quite often the transitions in both magnitudes take place simultaneously, 
except obviously when the QRMs reach the full spin-polarization, point from which on the 
isospin increases in one-unit jumps until the system is made up of only bonding states. 

The comparison of the isospin phases for QRMs with d = 4 and 6 nm is presented in Fig. 
E] for N = 32 and 40. Clearly, the highest values of I z appear for the smallest inter-ring 
distances, as expected from the single-particle levels shown in Fig. [31 corresponding to d = 4 
nm and N = 40, in which only a few antibonding orbitals are occupied for low values of 
B. Indeed, one can see from the bottom panels of Fig. [8] that for this inter-ring distance 
magnetic fields of about 5 T are enough to yield configurations with the maximum isospin 
value N/2, whereas for the QRMs with d = 6 nm such values of B still correspond to small 
J 2 's due to the large amount of filled antibonding states. 

We have also calculated the addition energies, defined by 

A 2 (iV) = E(N + 1)-2E(N) + E(N -1) , (3) 

E(N) being the total energy of the iV-electron system, for QRMs made of up to 14 electrons 
at different inter-ring distances, submitted to several magnetic fields, as a function of N. For 
the sake of comparison, we have also calculated A 2 (iV) for the corresponding single rings. 
The results for B = 0, 3, and 6 T are shown in Figs. [91TTH respectively, in which the bottom 
panels correspond to the single ring. 

From Fig. one can see that at zero magnetic field the single-QR addition spectrum 
presents the usual intense peaks at N = 2, 6 and 10 with zero total spin, and those at N = 4 
and 8 with S z = 1 satisfying Hund's rule. Similar results are found for the QRMs with 
d = 2 and 4 nm, indicating that such systems behave as a single ring owing to the strong 
quantum mechanical coupling corresponding to these inter-ring distances -notice that the 
spin values coincide for all the configurations but that with N = 13. This fact contrasts 
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where at d — 4 nm the 



with the results found for the vertical ring molecules of Ref. 
spectrum clearly reflected an intermediate coupling situation due to the filling of the first 
antibonding orbitals. As commented before, for the systems studied in the present paper 
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such states are only occupied for larger inter-ring separations (or iV's of the order of 30). 

The spectrum corresponding to d — 6 nm is shown in the top panel of the same figure. 
One can see that, although some of the marked peaks are preserved, in particular those at 
N = 2 and 8, the ones at N = 4, 6 no longer exist -notice that for 6 electrons the spectrum 
presents now a minimum and also that a new peak is found at iV = 5. This intricate structure 
can be understood from the corresponding single-particle energy levels. Indeed, it appears 
that the QRMs with N < 4 are made up of only bonding states, the first antibonding state 
being filled when N = 5. From N > 7 on, the QRMs have always occupied both B and AB 
orbitals but, however, the intermediate 6-electron configuration has again only symmetric 
states. This alternate behavior evidences that 6 nm is not a separation large enough for the 
QRMs to be in the weak coupling limit, but rather corresponds to an intermediate regime. 
Notice also that, from the results of Ref. [39!, in the weak coupling limit one would expect 



to find clearly marked peaks at the same N values as for the single ring multiplied by two 
-i.e. at N = 4, 12 and 20, indicating that the rings are so apart that behave as isolated 
entities. We have checked that for our QRMs to present such spectrum, we should consider 
inter-ring distances of about 10 nm. The different spin values for d = 6 nm as compared to 
those in the strong coupling regime can also be explained from the sp levels. For example, 
the 2S Z = 3 assignation of the QRM with iV = 5 is due to the above-mentioned filling of an 
antibonding -spin- up with / = 0- orbital replacing the spin-down \l\ = 1 state occupied for 
d = 2 and 4 nm. Analogously, the configuration with S z — 1 (instead of S z = 0) for iV = 10 
can also be explained from the sp levels: in the strong coupling limit, the QRM is formed 
by the spin-degenerated sp levels with I = 0, |1| and |2|, but this closed-shell configuration 
is prevented by the filling of the antisymmetric orbitals at d = 6 nm. Finally, the reverse 
situation occurs at N = 8, where the closing of the antibonding I = and |1| shells contrasts 
with the Hund's-rule configurations found for the strongly coupled molecules. 

Fig. [10] shows the addition energies corresponding to the situation in which a magnetic 
field of 3 T is applied to the rings. Like what is found at B = 0, the spectrum of the single 
system and of the molecules with d = 2 and 4 nm are rather similar -notice the different 
energy scales, the most remarkable difference being the salient minimum that appears for the 
single QR at iV = 5. For the above-mentioned inter-ring distances, peaks with S z = are 
found at iV = 2, 4, 8, 10 and 12, as well as a peak at N = 6 with S z = 2, although they are 
not as clearly marked as at B = 0. It turns out that, even in the presence of a magnetic field, 
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when the single-particle energy levels no longer display the ±Z— degeneration, the QRMs can 
adopt configurations that are somehow analogue to these characteristic of the situation at 
B = 0, namely the closed-shell ones and those fulfilling Hund's rule. Indeed, for e.g. d = 4 
nm and iV = 10, the ring molecule is made up of the spin-degenerate bonding states with 
I = 0—4 (instead of those with |/| = 0—2 of the B = case). Similarly, at iV = 6 the occupied 
orbitals are the spin- up and -down ones with I = 1 and 2, and the spin- up ones with / = 
and 3 (instead of the spin-degenerate states with \l\ = — 1 filled at zero magnetic field). 
For larger inter-ring separations, the occupancy of the first antibonding orbitals washes out 
these structures and the addition spectrum becomes flatter and irregular. One can notice 
also the different spin assignations between the single and the coupled systems, especially 
for the lowest-populated configurations. In particular, the single QRs with N < 5 turn out 
to be fully spin-polarized, which can be attributed to the combined effect of the magnetic 
field and a relatively strong exchange-correlation interaction characteristic of few-electron 
single quantum rings. The relatively higher spin values at d — 6 nm for N > 7 are due to 
the filling of the antibonding states. 

Finally, the addition energies for B = 6 T are shown in Fig. [TU It can be seen that 
in all cases the only clearly marked peak is the one at iV = 2, with the rest of the spectra 
being rather flat, following the trend observed at B = 3 T. Nevertheless, some weak peaks 
are still found and can be interpretated as in the previous cases, e.g. the one at N = 8 for 
d — 4 nm with 2S Z = 2: the system fills the spin- up and -down states with / = 2 — 4 and 
the spin-up ones with I = 1 and 5. One can also notice that the faint peak of the 4-electron 
configurations of both the single ring and the QRM with d = 2 nm becomes a minimum at 
larger inter-ring distances. Concerning the spin, the single QRs and the QRMs with d = 2 
and 4 nm turn out to be fully polarized for N < 7, 5 and 3, respectively, whereas the filling 
of the antibonding states favors the fully spin-polarization of molecules with the largest ring 
separation for all the considered electron numbers. 

IV. SUMMARY 

Within the local spin-density functional theory, we have addressed the ground state of 
quantum ring molecules containing up to 40 electrons, with different inter-ring distances, 
and submitted to perpendicular magnetic fields. In the strong coupling regime the energy 
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levels and the addition energies of the QRMs are similar to those of a single QR, although 
some differences are found due to the effect of the magnetic field, which has a tendency to 
wash out the clearly marked peaks characteristic of the B = case as well as to yield flatter 
addition spectra. However, even at B ^ 0, some peaks are still present and they can be 
interpretated as at zero magnetic field. 

When the ring separation is increased until the first antibonding orbitals are occupied, 
the addition spectra become irregular and the ring molecules are fully spin-polarized at 
relatively low magnetic fields. The filling of such states yields isospin oscillations as a 
function of B, increasing in one-unit jumps once the corresponding molecular configurations 
reach the maximum spin value. 

Despite the lack of experimental results to compare ours with, we believe that the ones 
herewith presented may be helpful in the analysis of future experiments on vertically coupled 
QRs concerning, e.g., the realization of single-electron transistor (SET) measurements, where 
the evolution of the chemical potential fi(N) with the magnetic field can be experimentally 
identified as the variation of the position of the current peaks as a function of the applied 
field, showing irregularities arising from phase transitions. 
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FIG. 1: Single-electron energies (meV) as a function of / for a N = 40 single ring. Up- 
ward (downward) triangles denote |(J.) spin states. The horizontal lines represent the Fermi levels. 
The value of B (T) is indicated in each panel. 
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FIG. 2: Same as Fig. [I] for a QRM with N = 40 electrons and d = 2 nm. Notice that due to the 
small separation between the rings only bonding states are occupied, the antbonding ones lying at 
much higher energies. 
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FIG. 3: Same as Fig. [2] for d = 4 nm. The increase of the inter-ring separation allows for 
the occupation of both bonding and antibonding states, indicated by solid and open triangles 
respectively. 
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FIG. 6: B values at which the selected QRMs become fully spin-polarized for the inter-ring 
distances d = 2,4, and 6 nm. The isospin of each configuration is also indicated. The lines have 
been drawn to guide the eye. 



22 




FIG. 7: Isospin (solid line, left scale) and spin (dotted line, right scale) values as a function of B 
for the QRMs with d = 6 nm and N = 8, 16 and 24. 
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FIG. 8: Isospin values as a function of B for the QRMs with d = 4 and 6 nm, and N = 32 and 40. 
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FIG. 10: Same as Fig. M for B = 3 T. 
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FIG. 11: Same as Fig. M for B = 6 T. 
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